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ABSTRACT 

Aims. We assess the ability of a solid ring to model a global perturbation induced by several thousands of main-belt asteroids. 
Methods. The ring is first studied in an analytical framework that provides an estimate of all the ring's parameters excepting mass. 
In the second part, numerically estimated perturbations on the Earth-Mars, Earth- Venus, and Earth-Mercury distances induced by 
various subsets of the main-belt population are compared with perturbations induced by a ring. To account for large uncertainties 
in the asteroid masses, we obtain results from Monte Carlo experiments based on asteroid masses randomly generated according to 
available data and the statistical asteroid model. 

Results. The radius of the ring is analytically estimated at 2.8 AU. A systematic comparison of the ring with subsets of the main belt 
shows that, after removing the 300 most perturbing asteroids, the total main-belt perturbation of the Earth-Mars distance reaches on 
average 246 m on the 1969 - 2010 time interval. A ring with appropriate mass is able to reduce this effect to 38 m. We show that, by 
removing from the main belt ~240 asteroids that are not necessarily the most perturbing ones, the corresponding total perturbation 
reaches on average 472 m, but the ring is able to reduce it down to a few meters, thus accounting for more than 99% of the total effect. 
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1. Introduction 

Asteroid perturbations are considered the major obstacle 
for achieving a satisfactory prediction accuracy in plane- 
tary ephemerides. The asteroid problem persists despite the 
ephemerides being fitted to a variety of highly accurate obser- 
vations available over almost 40 years. The most constraining 
observations are tracking data from spacecraft orbiting Mars, 
which over the past 8 years are available with metric preci- 
sion. Despite this accuracy, today's ephemerides are unable to 
extrapolate the Mars ranging data one year into the future with 



a precision better than 15 m (see Folkner et al. 2008 Fienga 



|et al. 2009 1. The obstacle to improving the situation is to cor- 
rectly account for a large number of similar asteroid effects with- 
out any accurate knowledge of the asteroid masses. In modern 
ephemerides, usually about 300 asteroids are modeled individu- 
ally and the rest of the main-belt is represented by a circular ring. 
Our objective in this study is to estimate the ability of a ring to 
model large numbers of asteroids. 

In section [2] we set a simple analytical framework and show 
that a ring is in fact a first-order representation of the main-belt 
perturbation. Section [3] describes a new implementation of the 
ring in the INPOP ephemeris (see Fie nga et al.||2009) . To eval- 
uate the ring model, we chose to compare the perturbation in- 
duced by the ring on the Earth-Mars distance with the perturba- 
tions induced by a test model containing thousands of asteroids. 
Section|4]describes the chosen test model that consists of 24635 
asteroids assigned each according to the available data with rea- 
sonable random mass values. In section|5 the effect of the ring is 
compared with the perturbation induced by the test model after 
removing various asteroid subsets. The comparison is repeated 
for different asteroid masses in Monte-Carlo-like experiments. 



We are thus able to estimate the number of asteroids that need to 
be individualized in the asteroid model so the ring is able to sat- 
isfactorily represent the remaining global perturbation. We also 
derive estimates for the amplitude of the global perturbation, to- 
gether with the corresponding mass of the ring and the ampli- 
tude of the residual perturbation unmodeled by the ring. The 
perturbations on the Earth-Mercury and Earth- Venus distances 
are considered at the end of the section. Some applications of 
the obtained results are discussed in section|6]. 

2. Analytical approach 

2.1. Averaging the main belt 

The perturbation induced on a planet by the main belt is the 
sum of the perturbations induced by the asteroids within the 
belt. Accounting in a model for each asteroid individually in- 
volves dealing with a very large number of unknown and highly 
correlated parameters. In consequence, it is not possible to im- 
plement all the individual main-belt asteroids in the model. We 
show that with some hypotheses the asteroid model can actually 
be reduced to the average effect of a single object. 

Let us consider a main-belt asteroid of mass m' on a fixed or- 
bit perturbing a planet of mass m. We denote the classical orbital 
elements with (a, A, e, m, I, £2). In this whole section we use the 
convention of marking the variables related to the asteroid with 
a prime and the variables related to the planet without a prime. 
We assume that in the main-belt the considered asteroid is not 
alone on its orbit, but many other objects with similar masses 
are spread uniformly in terms of mean longitude on similar or- 
bits. It is then possible to average the perturbation of the consid- 
ered asteroid over the mean longitude A' without losing any part 
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Table 1. Parameters of an averaged orbit representing the main- 
belt perturbation as determined from the perturbation on the in- 
ner planets. 



planet 


a [AU] 


e' 


m' n 


/' [°] 


fi' [°] 


Mars 


2.71 


0.04 


-1.30 


0.96 


-60.34 


Earth 


2.74 


0.04 


-4.88 


0.89 


-50.31 


Venus 


2.75 


0.04 


-4.03 


0.89 


-51.08 


Mercury 


2.75 


0.04 


-5.71 


0.88 


-50.67 
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To obtain the perturbation induced by the entire main belt, the 
above equations have to be summed over all the asteroids. For A 
we can write 



= f j X i = N m (n^(abZ-b%) 



of the asteroid's contribution to the main-belt perturbation. For 
the sake of simplicity, we also average the perturbation over the 
mean longitude of the planet. We are thus left only with a secular 
effect. 

Laskar & Robutel ( 1995 ) provide the secular Hamiltonian of 



the three-body problem expanded in eccentricity and inclination. 
In the case where the semi major axis of the asteroid is always 
greater than the semi major axis of the planet, the Hamiltonian 
can be rewritten up to the second degree as 



F(K,A,H,h,Z,0 
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where b\j 2 , by 2 , by 2 are Laplace coefficients depending on the 
semi-major axis ratio a = a/a' < 1. The canonical variables 
(A, A, H, h, Z, are obtained by a linear transformation from the 
variables of Delaunay: 

A = m -\JjIa 



H = mv£a(l - Vl - e 2 ) 
Z = m Jfia(l - e 2 )( 1 - cos /) 



h = -m, £ = -£X 

We can calculate the perturbation induced on the planet by the 
asteroid from Hamilton's equations. By only keeping the lowest 
order terms in eccentricity and inclination, we get 

a = 



A = na2 'i^( ab f/ 2 ~ b v 2 ) 
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where () represent the average over all the asteroids, and N m 
denotes the total number of modeled objects. If we assume that 
there is no correlation between mass and orbital elements we 
obtain 



X tot = N m { m ')l^(abf/ 2 -b^) 



Thus A lo , can be obtained by averaging the A term with m' = 1 
in equations ([TJ and multiplying the result by the total mass of 
all the main-belt asteroids denoted here by m' tot . The orbits of 
the asteroids are determined fairly well. We use asteroid orbits 
available in the ASTORBQdatabase to calculate 



1 N - 2 

N m jrf m e ^ 3 / 2 



numerically for all the inner planets. We consider in the cata- 
log only those asteroids with absolute magnitude below 14 and 
semi-major axis below 3.5 AU. Thus the estimation is based on 
a total of 24634 orbits. There are actually 24635 asteroids that 
satisfy these criteria. We eliminated 433 Eros from the selection 
because its semi-major axis is lower than the semi-major axis of 
Mars and thus equations ([T]i do not apply. 

The method described to calculate A lo ,/m' tot is applied also 
to the other orbital parameters. We thus obtain the numerical 
values of e m /m' tot , vy tot jm' m , i tot /m' tot and Q, tot /m' tot . For each 
inner planet, we can find an average equivalent orbit that will 
perfectly fit all the calculated total perturbations. Because A is 
a monotonously decreasing function of a, the value of k tot /m' tm 
determines the semi-major axis of the equivalent orbit unam- 
biguously. Once the semi-major axis is determined, the calcu- 
lated values of e tot /m' tot , vj t otlm' tot determine the eccentricity and 
perihelion. Similarly i, or /m' ton Q. ro ,/m' tot determine the inclina- 
tion and node. Table [T] summarizes the orbital parameters of 
the equivalent orbit obtained on all the inner planets by aver- 
aging the 24634 asteroid orbits. The chosen reference frame is a 
nominal invariable plane defined by an inclination of 23°00'32" 
and a node at 3°51'9" with respect to the International Celestial 
Reference Frame. Equations ([T} contain only first-order terms in 
eccentricities and inclinations. Because of this approximation, 
the parameters in Table [T] vary from one planet to another. 

The effects of an orbit averaged over its mean longitude are 
entirely equivalent to the effects induced by a solid ring. This 
equivalence has been known since the works of Gauss (see for 



The asteroid orbits are calculated for the 27 October 2007. 
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example |Hill| |1882| l and is valid in both the secular and non- 
secular cases. For practical purposes, we refer in the following 
more often to a ring than to the effect of an averaged orbit, but 
the reader has to bear in mind that both are equivalent. 

Although we observe non-zero eccentricity values in Table 
[T[ in this work we study the modeling of the main belt with a 
circular ring. The non-zero eccentricities are a consequence of 
the tendency of asteroid perihelia to accumulate close to the lon- 
gitude of the perihelium of Jupiter. The precession rate of an 
asteroid perihelium is in fact at its lowest when close to the per- 
ihelium of Jupiter. The asteroid perihelia thus spend more time 
in this region, which explains their apparent accumulation (see 
|Murray & Dermott 1999| l. For the remainder of this paper we 
fix the radius of the ring at 2.8 AU. This value is within 0. 1 AU 
of all the radii in Table [T] and it is also the value of the radius 
of the ring adopted in INPOP06 ( |Fienga et al] [2008| . The ini- 
tial inclination of the ring with respect to the invariable plane is 
chosen at zero. This fixes all the ring's geometrical parameters. 
The only parameter of the main-belt model that is left undeter- 
mined is the mass of the ring or equivalently the total mass of 
the modeled asteroids. 



2.2. Perturbation of the orbital elements 

We use equations (JTJ to obtain analytical expressions of the sec- 
ular perturbation induced by a ring on a planet. We stress that, 
although the analytical results are presented in the secular case, 
it is merely for convenience and the objective in this study is to 
assess the ring model in the general non-secular case. By setting 
e' and /' at zero one obtains 

a — e — I — 



Table 2. Secular effect of a ring, with INPOP06 parameters, on 
the orbital elements of the Earth and Mars. 
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0.94 
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i 
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rad yr~' 


a. 


-0.94 


-2.65 


rad yr~' 



The table shows that the secular effects on the eccentricities and 
inclinations are clearly negligible. 

2.3. Perturbation of the Earth-Mars distance 

In this work we focus on the Earth-Mars mutual distance, which 
today is the most accurately observed parameter. Indeed Mars 
has been a target for many missions. These provide highly accu- 
rate ranging data with accuracies varying from roughly 20 m for 
the Viking data to accuracies of about 1 m for more recent mis- 
sions like Mars Global Surveyor, Mars Orbiter, or the ongoing 
Mars Express mission. We limit our study of the asteroid pertur- 
bations to the interval between years 1969 and 2010. This corre- 
sponds roughly to the interval spanned by the available ranging 
data. 

Let us denote with AD any perturbation induced on the 
Earth-Mars distance. We have 



,m 



TIT 



o 



2 m tot , (1) 

na - — b\L 
Atria 1 



i m tot h m 

4m 3 / r 



(2) 



The ring induces linear drifts in the mean, perihelion, and node 
longitudes of the planet. There is no secular effect on the semi- 
major axis because, in the secular case, the Hamiltonian is in- 
dependent of A. The presence of the ring may change the mean 
semi-major axis by a small constant value Aao- This small value 
would translate as a shift in the mean motion Ann (= -3nAaa/2a) 
that should be accounted for in A. Equations (Kb predict no per- 
turbation of eccentricity or inclination for the second-degree 
Hamiltonian. To estimate the effect on these parameters, we need 
to consider at least a Hamiltonian of degree 4. With Laskar & 
|Robutel| ( p"995| > we then obtain 



= no 2 '^-I 2 e sin (2m - 2Q) {-—abfl 

m Q \ 32 3/z 
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In Table [2] we use equations (j2j) and Q to quantitatively es- 
timate the effect induced on Earth and Mars by a ring with radius 
at 2.8 AU and a mass of 0.34 x 1O~ 1O M . These are parameters 
of the asteroid model in INPOP06 where the ring represents all 
the main-belt asteroids except for a selection of 300 individuals. 



AD = D - D , 

where D and Do are the Earth-Mars distances in the perturbed 
and unperturbed cases respectively, and D and Do each depend 
on the perturbed and unperturbed orbital elements of Earth and 
Mars. By writing a first-order Taylor expansion of D in terms of 
the perturbations of the individual orbital elements, we obtain 



AD 
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(4) 



The indexes E and M refer to the Earth and Mars. In the expan- 
sion, (Aa, AA..., AQ.) represents the perturbations of the planetary 
orbital elements. We have, for example, for the perturbation in 
mean longitude 

AA - A - A 

where A and An are the evolutions of the mean longitude in the 
perturbed and unperturbed cases, respectively. In an analogous 
way, we can define the perturbations of the other orbital ele- 
ments. 

We calculated estimates of the partial derivatives by differen- 
tiating a second-degree eccentricity and inclination expansion of 
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Table 3. The Earth-Mars distance dependence on the orbital el- 
ements. 
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the Earth-Mars distance. The resulting estimates are functions 
of approximately constant amplitude oscillating with frequen- 
cies close to the Earth-Mars synodic frequency. Table [3] lists the 
amplitudes of all the partial derivatives. The significant differ- 
ences among the amplitudes of the partial derivatives stem from 
the various orders in eccentricity and inclination. 

The general expression Q can be used to estimate the effect 
induced on the Earth-Mars distance by the ring with INPOP06 
parameters. With Table [3] and the values of the secular pertur- 
bations calculated in Table |2j we find that the perturbation will 
reach approximately 5 m over one year, mainly because of the 
drift in the mean longitude of Mars. 



3. Numerical integrations 

3.1. The implementation of the ring in INPOP 

In Krasinsky et al.| p002 ) the ring was numerically implemented 
as a perturbing acceleration in the ecliptic plane of each planet, 
and this implementation has been adopted in INPOP06 as well. 
Fienga et al. ( 2009) show that in long-term integrations the ring 



causes a slight drift in the Solar System barycenter of approxi- 
mately 10 m over a century. Indeed with the ring modeled as an 
exterior force, the system's total linear momentum is not con- 
served. The existence of this drift motivated a more realistic im- 
plementation of the ring adopted since in INPOP08. 

The ring is treated as a solid rotating body, which fully inter- 
acts with the planets. Although its initial orientation is taken to 
be parallel to the system's invariable plane, the orientation of the 
ring is an integrated parameter that evolves with time under the 
influence of the planetary perturbations. The angular momentum 
of the ring is constant in amplitude and determined by the radius 
of the ring and Kepler's law. The linear and angular momenta of 
the system are thus conserved, which eliminates the barycenter 
drift occurring in the previous implementation. Because the av- 
eraged orbit is bound to the Sun by gravitation, in INPOP we fix 
the center of the ring to the barycenter of the Sun. We note that a 
free floating ring centered on the Sun or on the barycenter of the 
Solar System is actually unstable and would gradually drift away 
from its initial position. The expression for the force exerted by 
a ring on a particle is given in Appendix [A] 

3.2. Comparison with the analytical expressions 

We used INPOP to numerically estimate the effects of a ring with 
INPOP06 parameters on the Earth and Mars. These effects can 
be isolated from the perturbations induced by other objects in the 
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Fig. 1. Perturbation induced on the orbital elements of Earth (in 
gray) and Mars by a ring with INPOP06 parameters. 



INPOP model by comparing a Solar System integration with the 
ring and a reference integration in which the ring is absent. The 
difference between the evolutions of orbital elements in both in- 
tegrations provides the perturbation induced exclusively by the 
ring. Results obtained with this method on the 1969-2010 inter- 
val are shown in figure [TJ Because integrations both with and 
without the ring start from identical initial conditions at J2000, 
all the perturbations are at zero for this date. A linear regression 
provides values for the numerically observed secular drifts : 

= -1.15 x 10-" radyr 
= 0.95 x 10" radyr 1 



..-1 



X M = -3.11 x 10"" radyr 
m M = 2.65 x 10~ n radyr* 1 



-0.94 x lO" 11 rad yr 1 , Q. M = -2.69 x 10" 11 rad yr 



With the exception of the mean longitude of Mars, the observed 
drifts agree well with the analytical predictions of Table [2] For 
Mars, the discordance with the predicted drift in mean longi- 
tude is due to the non zero mean value of the perturbation of the 
semi-major axis of Mars. A constant shift in the semi-major axis 
corresponds to a shift in mean motion A«o- We can estimate A«o 
directly from the mean value of the semi-major axis in figure 
[T] Adding the estimated A«o to the drift in the mean longitude 
of Mars calculated in Table |2] leads to a corrected estimation 
Am = -3.05 x 10~ n rad yr -1 . The numerically observed drift 
in the inclination of Mars is 7 M = 0.0024 x 10~" rad yr 1 . The 
interval 1969-2010 is too short to correctly estimate the secular 
drift in the inclination of Earth. Similarly, for both planets it is 
impossible to estimate the secular drifts in eccentricities. 

We estimate the perturbation induced by the ring on the 
Earth-Mars distance by comparing two integrations with and 
without the ring. Figure|2]shows that with INPOP06 parameters 
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AD [ m ] 



Table 4. Albedo and corresponding uncertainty for objects with 
known taxonomy or belonging to a dynamical family. 




Fig. 2. Perturbation of the Earth-Mars distance induced by a ring 
with INPOP06 parameters. 



the effect of the ring reaches approximately 150 m. This pertur- 
bation results from secular drifts in the various orbital elements, 
thus the 150 m reached over 31 years between 1969 and 2000 
can be translated as approximately 5 m per year. This is in good 



agreement with the prediction made at the end of section 2.3 



4. A test model of the main belt 

From the analytical standpoint the ring appears as a first-order 
representation of the main belt. Nevertheless, it is rather diffi- 
cult to evaluate its effective capacity to model a multitude of 
objects. One possible solution is to test the ring against a test 
model containing large numbers of individual asteroids assigned 
with particular masses. The ring can be tested against many dif- 
ferent asteroid models each built with a different set of masses. 
Such Monte Carlo (MC) experiments can provide an estimate of 
the number of asteroids that have to be modeled individually in 
ephemerides. If the asteroid masses are individually inaccurate 
but globally reasonable, the MC experiments can also estimate 
the ring's mass. 

4.1. Asteroid selection 

We select for our test model all the asteroids in the ASTORB 
catalog that have an absolute magnitude brighter than 14 and 
are situated in terms of semi-major axis within 3.5 AU. These 
criteria are the same as in section |2.1| The absolute magnitude 
limit of 14 corresponds to the estimated completeness limit of 
the main-belt and NEO populations as reported by Jedicke et al. 
(2002 ). The value also leads to a relatively large but still reason- 
able number of objects to work with, 24635 asteroids in total. 
Absolute magnitude can be converted to diameter with the fol- 
lowing formula (Bo well et al.||1989] l: 



(km) 



1329 , 



■10 



-0.2H 



where D is the diameter in kilometers, H the absolute magnitude, 
and p the albedo of the asteroid. With a value for a minimum pos- 
sible albedo within the asteroid population, the formula sets an 
upper bound on the diameter of asteroids that are not included in 
our model. [Tedesco et aIT| ( |2005( l estimate the mean albedo of the 
lowest albedo class at approximately 0.05. With this estimate, 
the maximum diameter of an asteroid not included in our test 
model is 10 km. 
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Dora 


0.0603 


0.0160 


Eos 
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0.0150 
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Taxonomy data " 




Class 


P 


Ap 


Low albedo (C,G,B,F,P,T 


D) 0.0545 


0.0345 


Intermediate albedo (M) 


0.1005 


0.0155 


Moderate albedo (S,Q) 


0.2335 


0.1215 


High albedo (E,V,R) 


0.4305 


0.0955 



Values reproduced from Tedesco et al. (2005 1, Tables 1 and 7 



4.2. Asteroid masses 

Today the total number of accurately measured asteroid masses 
amounts to only a few tens. Besides the case of binary objects, 
asteroid mass determinations are in general susceptible to sys- 
tematic errors, which are hard to estimate. To assign all the se- 
lected asteroids with at least reasonable masses, we devise a sim- 
ple algorithm inspired by the statistical asteroid model ( |Tedesco 
et al. 2005) >. The algorithm processes family membership, tax- 
onomy, and SIMPS survey dat£0 and assigns accordingly each 
asteroid with an albedo p and corresponding uncertainty Ap. 
Among the 24635 asteroids, more than two thousand asteroids 
have values of p and Ap directly available from SIMPS, objects 
with family or taxonomy data have their p and Ap attributed with 
data obtained in the statistical asteroid model and reproduced 
in Table [4] In cases where information overlaps, SIMPS is pre- 
ferred to taxonomies and families. Taxonomy data is preferred 
to families, which helps to eliminate interlopers in the family 
data. In general at least some information is available for 10% 
of the selected objects. For the remaining majority of asteroids, 
they are randomly assigned an albedo class with probabilities 
adopted from Tedesco et al. (2005 ) (56% low albedo, 7% in- 
termediate albedo, 34% moderate albedo, and 3% high albedo). 
Table |4]defines p and Ap for each albedo class. 

Asteroid diameters are calculated from corresponding p to 
which a random error within ±Ap is added. For SIMPS data, 
the algorithm ignores the formal Ap and applies instead a more 
realistic 10% uncertainty directly on the diameter. As potential 
systematic errors in absolute magnitudes are reported by various 
authors ( Juric et al. 2002 ), we account for them with a ±0.5 ran- 
dom uncertainty added to absolute magnitude. With the lower 



The databases are maintained in NASA's PDS Asteroid Archive. 
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part of Table|4j we can use each asteroid's mean albedo p to as- 
sign the asteroid with a density class : C (low albedo), S (mod- 
erate albedo), and M (intermediate and high albedos merged). 
We stress that these density classes are attributed according to 
albedo. Therefore, for some objects, the density classes may not 
coincide with taxonomy data. Bulk porosity is expected to vary 
among the asteroids (see |Britt et al.| [2002 1, so we adopt the fol- 



[m] 



lowing intervals for the class densities : [0.5,2.5] (C), [1.6,3.8] 
(S) and [1,5] (M). The density of an asteroid is chosen randomly 
within the corresponding class interval and together with the pre- 
viously calculated diameter provides a mass. The masses of six 
asteroids are kept constant and equal to their published values. 
We fix 4.756 x 10- 10 M o , 1 .025 x 1Q- 10 Mq, 1 .348 x 10" 10 M Q for 
1 Ceres, 2 Palla s, 4 Vesta (|Fienga et ak] |2008|l, 0.45 x 10" 10 M G 
for 10 H vgiea (jChesley et al.l |20 = 0y 0.03 x 10- 10 M o for 22 
Kalliope (Mer line et al.| |1999), and 0.037 x 10~ 10 M Q for 45 
Eugenia ( |Margot & Brown||2005) . 

We used the algorithm to generate for each asteroid a set of 
100 random masses. A standard mass set is generated without 
any random choices. In this standard mass set, Ap is put to zero 
and densities are maintained at INPOP06 values : 1.56 (C), 2.18 
(S), and 4.26 (M). Objects without any available data are auto- 
matically considered as belonging to the C taxonomy class in the 
standard mass set. 



4.3. Individual perturbations and the global effect 

As in section [3] the perturbation of the Earth-Mars distance is 
denoted with AD. To evaluate the perturbations AD, (1 < i < 
24635) induced by each individual asteroid of the test model, 
we performed extensive integrations with INPOP. Each AD, was 
obtained by comparing on the 1969-2010 interval a Solar System 
integration with the particular asteroid and a reference integra- 
tion in which the asteroid is absent. 

For a given set of asteroid masses, we can rank the aster- 
oids according to the decreasing amplitude of their individual 
perturbations. Here and in the following, the amplitude of a per- 
turbation is estimated by the maximum of |AD| reached on the 
1969 - 2010 interval. Each AD, is proportional to the mass M, of 
the perturbing asteroid, in consequence, we have 

AD,- = d 4*M, (5) 

Equation |5]l can be rewritten as 
dD 

AD,- = —Mi 
dMi 

where D is the Earth-Mars distance. This is the analog of equa- 
tion but instead of considering the perturbation of the Earth- 
Mars distance as depending on the perturbations of the planetary 
orbits, we consider it as depending on the mass of the perturbing 
asteroid. The first-order Taylor expansion of D in terms of all the 
asteroid masses leads to 



D(M 1 ,...,M 246 35) = O(0, 



dD dD 
, 0) + - — Mi + ... + 

<9Mi <9M 24 635 



M24635- 



We can thus approximate to first-order the perturbation caused 
by a particular set of asteroid masses by the sum of the already 
calculated individual asteroid contributions: 



AD = ADi + ... + AD24635- 



(6) 



Although initially the AD, are obtained with INPOP for the stan- 
dard mass set defined in section |4~2] the proportionality relation 




Fig. 3. Difference between the AD obtained from expansion (Jf^jl 
and the AD obtained from a simultaneous INPOP integration of 
the asteroids. 

Q provides the individual perturbations for any set of aster- 
oid masses. In the next section we will use the expansion (|6]l 
to study how the resemblance between the ring, and the test 
model evolves with the number of asteroids removed from the 
test model and between different mass sets. Indeed, INPOP in- 
tegrations of the Solar System with thousands of asteroids take 
several hours, and it is impossible to explore the parameter space 
by reintegrating the test model each time with new parameters. 
We refer to the perturbation induced by the test model after re- 
moving a selection of individuals as to the global perturbation. 

To test the development (|6j, we can compare the global per- 
turbation obtained from a simultaneous INPOP integration of all 
the asteroids, but the N most perturbing ones with the same per- 
turbation obtained from the development. We choose N at 300 
here because it corresponds to the number of asteroids usually 
considered individually in modern ephemerides. Figure[3]shows 
the difference between the two perturbations for the standard 
mass set. This difference is on the order of 1 m, which is less than 
1% of the perturbation's amplitude. In consequence, we consider 
the development (|6]) as satisfactory. It should be noted, that in the 
simultaneous integration, the mutual perturbations of the aster- 
oids were not taken into account. 



5. Testing the capacity of the ring to model large 
numbers of asteroids 

5. 1 . Selection based on amplitude 

We denote with AD g i ob (A0 the global perturbation induced on 
the Earth-Mars distance by the test model after removing from 
the test model the N most perturbing asteroids. Similarly we de- 
note with AD r i ng the perturbation induced by a ring. To evaluate 
the capacity of the ring to represent the global perturbation, we 



fit for different values of N the amplitude of AD r i ng so as to min- 
imize 

|AD ring - AD glob (A0| . 

Because AD r i ng is proportional to the mass of the ring, fitting 
the amplitude is actually equivalent to fitting the mass of the 
ring. Let us denote by R(N) the percentage of the amplitude of 
the global perturbation left after fitting the ring. We thus have 



R(N) 



max |AD ring - AD glob (A0| 
max|AD glob .(A0| 
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Fig. 4. Evolution of R(N) for the standard set of masses (in gray) 
and an average over 100 different sets. 

Figure |4] shows the evolution of R(N) for the standard set of 
masses, as well as an average of R(N) over the 100 different mass 
sets defined in section 4.2 For values of N greater than 200, the 



ability of the ring to model a global effect has reached its max- 
imum and remains constant. At its best, the ring is thus able to 
represent more than 80% of the total perturbation amplitude. 

We show in figure [5] the effect induced on the planetary or- 
bital elements by the test model after removing from the test 
model the 300 most perturbing asteroids (for the standard mass 
set). The observed effects are indeed rather smooth and similar 
to the drifts induced by a ring. For comparison, figure [6] shows 
the effect induced on the planetary orbital elements by all the 
asteroids of the test model. 

It is possible to calculate the average mass of the ring after 
removing the 300 most perturbing asteroids by minimizing 

lADrfng - AD glob (300)|. 

With the 100 different sets of asteroid masses, the mass of the 
ring is estimated at M r - mg = 0.6 + 0.2 x 10 -10 M o . In a similar 
way, we can estimate the average maximum reached by AD g \ \, 
at 246 m and the average maximum reached by the residuals af- 
ter fitting the ring at 38 m. The obtained value of M r j ng is approx- 
imately twice as high as the mass of the ring fitted in INPOP06. 
Contrary to INPOP06 where parameters are fitted to observa- 
tions, here the value of M r - mg is obtained solely from the test 
model. Figure [7] shows the variations with N of the amplitude of 
the global perturbation and the amplitude of the residuals. The 
evolution of the mass of the ring is not shown because it follows 
proportionally the amplitude of the global perturbation. 

The calculated residuals for N = 300 correspond to approx- 
imately 10 m over the 2000-2010 time interval. This is an order 
of magnitude above the residuals obtained today for the most ac- 
curate Mars ranging data. The limiting factor of the ring model 
is the inability to reproduce the quadratic evolution in the mean 
longitude of Mars (see figure [5]). This quadratic evolution is in 
fact a consequence of the linear drift of the semi-major axis of 
Mars that persists in the test model for all values of N. 

5.2. The selection as a mixed integer quadratic problem 

The progressive removal of asteroids from the test model can 
be understood as a selection of individual asteroids that should 
be modeled individually in an ephemeris. The selection scheme 
based on the amplitude of the individual asteroid perturbations 
on the Earth-Mars distance is not optimal. Indeed for each state 
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Fig. 5. Perturbation induced on the orbital elements of Earth (in 
gray) and Mars by the test model after removing from the test 
model the 300 most important perturbers. 



N, it is possible to slightly modify the set of the removed aster- 
oids in order to eliminate the linear drift in the semi-major axis of 
Mars responsible for the quadratic evolution of the mean longi- 
tude in figure[5] The modification consists of removing from the 
test model a few additional asteroids that induce a positive slope 
in the perturbation of the semi-major axis of Mars and adding the 
same number of already removed objects with a negative slope. 
Such changes improve the average residuals from the previous 
38 m to 20 m. 

A more systematic approach is to use a combinatorial opti- 
mization algorithm to select among the N most perturbing aster- 
oids those that should be removed from the test model in order 
to maximize the modeling capacity of the ring. The problem can 
be stated formally for a particular mass set and a given N as the 
search for N+l parameters a, that minimize 

|a,ADi + ... + a N AD N + a N+ iAD dng + AD glob (A0| 

with the constraint that a N+1 is real positive and all the other a, 
are binary (equal to 1 or 0). The problem falls into the category 
of mixed integer quadratic problems (MIQP). In this particular 
case, a direct search for a solution involves testing 2 N combina- 
tions and becomes very difficult for any useful value of N. The 
exact solution can nevertheless be found with various methods, 
Gueye & Michelon ( 2009 \ provide a short literature review on 



the subject. One method of solving the MIQP is to linearize it 
and solve the linear formulation with the simplex algorithm pro- 
vided for example in the GNU Linear Programming Kit (GLPK). 
A classical linearization is Glover's linearization (see Gueye & 
Michelon, 2009), which involves N additional variables and 2N 
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Fig. 6. Perturbation induced on the orbital elements of Earth (in 
gray) and Mars by the entire test model. 
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Fig. 7. Evolution of the amplitude of AD g \ \, and the amplitude 
of the residuals after fitting a ring for the standard set of masses 
(in gray) and corresponding averages computed over 100 ran- 
dom mass sets. 



additional linear constraints. Unfortunately its implementation 
and the subsequent solution with GLPK does not lead to a solu- 
tion within any reasonable time. Experiments with the commer- 
cial optimization package CPLEXn on the other hand, showed 
an ability to solve the quadratic problem within seconds for N 
below 120 and within minutes for N below 200. For higher val- 



1 The software solves a hierarchy of linear subproblems in a branch- 
and-bound approach, see www.ilog.com/products/cplex for more de- 
tails on the package 



Fig. 8. Evolution of R(N) with the MIQP selection for the stan- 
dard set of masses (in gray) and an average over 100 different 
sets. The dashed line represents the average R(N) obtained with 
the selection based on amplitude. 
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Fig. 9. Evolution with the MIQP selection of the amplitude of 
ADgi b* and the amplitude of the residuals after fitting a ring 
for the standard set of masses (in gray) and corresponding av- 
erages computed over 100 random mass sets. The dashed lines 
represent the average amplitudes and residuals obtained with se- 
lection based on amplitude. 



ues, the time taken to solve the problem becomes too large to be 
of practical use in our MC experiments. 

Figure [8] gives the analog of R(N) obtained with a new 
selection scheme using the MIQP formulation, we denote the 
global effect corresponding to this new scheme with AD i b*(N). 
Because the calculations are relatively time consuming, only val- 
ues of N below 500 are considered. For each N the algorithm 
uses CPLEX to select among the N most perturbing asteroids, 
those that should be removed from the test model to obtain an 
optimal fit with a ring. For N greater than 200, only the 200 least 
perturbing asteroids among the N most perturbing ones are con- 
sidered by CPLEX. The remaining (N-200) are removed from 
the test model automatically. In figure [8] we observe that the 
ring's modeling capacity is greatly improved when compared 
with the selection scheme based on amplitude. Not only does 
R(N) drop almost to zero, but the maximum modeling capacity 
is also reached earlier. Figure [9] shows the evolution of the av- 
erage maximum reached by the global perturbation and of the 
corresponding residuals after fitting the ring. Although the am- 
plitude of the global effect is almost twice as large as in figure 
[7] the residuals are greatly improved. The performance reached 
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Fig. 10. Perturbation induced on the orbital elements of Earth 
(in gray) and Mars by the test model after removing from the 
test model at most 300 asteroids with the MIQP selection. 



after removing the 300 most perturbing asteroids from the test 
model is obtained with the MIQP approach for N = 50. The 
average residuals after fitting the ring are 4 m for ^ = 300. 
This corresponds to approximately 1.3 m over 10 years. The 
unmodeled part of the global perturbation is thus below the 
one sigma residuals of the MGS/MO and MEX data obtained 
in INPOP08. The average maximum reached by the global ef- 
fect is 472 m, which corresponds in terms of the ring's mass to 
Mnng = 1.1 +0.2x 10- ll) M Q . 



We show in figure 10 the perturbation induced on the or- 



bital elements of Earth and Mars by the global effect obtained 
with the MIQP selection (N = 300) and the standard set of as- 
teroid masses. The quadratic evolution of the mean longitude of 
Mars from figure 5] is straightened up. This explains why the 
corresponding global effect on the Earth-Mars distance has in- 
creased in amplitude. The MIQP approach in general improves 
the resemblance between the effect on the orbital elements in- 
duced by the global perturbation and the ring. A comparison be- 
tween figures [T] [5] and 10 shows that beside improvements in 
the semi-major axis and mean longitude perturbations, the new 
selection also improves the eccentricity and perihelion perturba- 
tions. There are still some discrepancies, most importantly in the 
perturbations induced on the inclinations and nodes as well as on 
the perihelium of the Earth; however, these discrepancies are not 
surprising as the asteroid selection is based only on the Earth- 
Mars distance, which according to Table [3] is not very sensitive 
to these parameters. 




Fig. 11. Average residuals computed with the MIQP selection 
over 100 random mass sets for the Earth-Mars (a and a'), Earth- 
Venus (b and b'), and Earth-Mercury distances (c and c'). The 
continuous lines represent residuals obtained from a simultane- 
ous fit of the ring to the AD„] b* of the three planets, and the 
dashed lines represent the residuals obtained by fitting only one 
A^ulob* at a time as in section |5. 2 



5.3. Accounting for all the inner planets 

Today ephemerides are being fitted to accurate Venus rang- 



ing observations from the ongoing VEX mission (see Fienga 
|et al.| |2009[ ). According to |Ashby et al. ( 2007| l, Mercury ranging 
could become available within a few years with missions like 
Messenger or BepiColombo. To evaluate the ability of the ring 
to model the global effect on all the inner planets, we repeated 
the MC experiment made with the MIQP selection scheme in 
Instead of fitting only the Earth-Mars distance, we 



5.2 



section 

simultaneously fit the effects on distances to Mercury, Venus, 
and Mars. Figure 11 shows the evolution of residuals for all the 
planets in the simultaneous fit. For comparison the figure shows 
the residuals that can be obtained by fitting the distance to each 
planet separately. For N = 300, the ring is able to represent the 
global effect simultaneously on all the inner planets with an av- 
erage accuracy better than 1.6 m over a 10 years time interval. 
Fitting all the planets together for N = 300 leads to the same 
amplitude of the global perturbation as in section 5.2 and hence 
to the same estimate of the ring's mass. When considering sep- 



arately Mercury or Venus, figure 1 1 shows that residuals well 



below 1 m per year can be reached by removing merely 1 Ceres, 
2 Pallas and 4 Vesta from the main belt. 



6. Applications 

6. 1 . Model selection 

To model the global effect correctly with respect to the accuracy 

is 



of available data, a reasonable choice according to figure 1 1 



to account for fewer than 300 individual objects. By examin- 
ing asteroids that were removed from the test model during the 
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Table 5. Asteroids with an effect on the Earth-Mars distance greater than 100 m. The last column corresponds to the probability of 



being removed from the test model during the Monte Carlo experiments (see section 6. 1 



Asteroid name 


Diameter [km] 


Density class 


a 


Earth-Mars [m] 


Earth- Venus [m] 


Earth-Mercury [m] 


prob. [%] 


4 Vesta 


468.30 


c 


11455.34 


164.42 


323.36 


100 


1 Ceres 


848.40 


c 


7829.72 


2656.87 


935.08 


100 


2 Pallas 


498.06 


c 


6612.41 


977.16 


379.08 


100 


324 Bamberga 


229.44 


p 


1820.88 


17.04 


19.46 


100 


10 Hygiea 


407.12 


p 


806.70 


132.18 


61.95 


100 


6 Hebe 


185.18 


c 


691.70 


54.66 


21.15 


100 


532 Herculina 


222.38 


c 


660.82 


24.04 


6.53 


100 


19 Fortuna 


213.47 


p 


572.41 


18.53 


9.37 


100 


51 Nemausa 


147.86 


1V1 


504.89 


33.92 


12.06 


100 


7 Iris 


199.84 


c 


489.73 


37.20 


12.85 


100 


9 Metis 


152.70 


c 


364.56 


4.82 


1.87 


100 


20 Massalia 


145.50 


c 


285.84 


20.90 


7.31 


100 


139 Juewa 


156.60 


p 


281.71 


5.07 


2.78 


100 


31 Euphrosyne 


255.90 


p 


279.07 


18.07 


7.00 


100 


16 Psyche 


253.16 


c 


264.08 


38.66 


13.92 


100 


52 Europa 


302.50 


p 


252.60 


43.76 


15.93 


100 


24 Themis 


218.44 


p 


234.96 


2.07 


1.25 


100 


63 Ausonia 


103.14 


c 


234.40 


4.15 


1.44 


100 


747 Winchester 


171.72 


p 


222.34 


27.55 


13.06 


100 


60 Echo 


60.20 


c 


165.80 


0.15 


0.15 


100 


78 Diana 


120.60 


p 


157.41 


1.49 


0.41 


100 


354 Eleonora 


155.16 


c 


152.73 


0.37 


1.88 


100 


41 Daphne 


174.00 


p 


149.44 


6.11 


4.25 


100 


3 Juno 


233.92 


c 


145.29 


15.88 


5.58 


100 


654 Zelinda 


127.40 


p 


142.71 


21.79 


13.53 


100 


5 Astraea 


119.06 


c 


139.25 


10.29 


3.93 


100 


18 Melpomene 


140.56 


c 


130.84 


18.96 


14.11 


100 


128 Nemesis 


188.16 


p 


128.69 


17.32 


7.71 


100 


192 Nausikaa 


103.26 


c 


128.46 


0.67 


2.26 


100 


11 Parthenope 


153.34 


c 


115.67 


5.73 


3.86 


100 


106 Dione 


146.60 


M 


114.27 


4.99 


1.58 


100 


105 Artemis 


119.08 


c 


113.20 


1.33 


0.59 


100 


23 Thalia 


107.54 


S 


1 12.92 


1.05 


0.56 


100 


372 Palma 


188.62 


C 


112.84 


3.79 


1.58 


100 


29 Amphitrite 


212.22 


s 


106.71 


43.13 


13.78 


100 


419 Aurelia 


129.00 


c 


106.21 


3.20 


1.23 


100 


15 Eunomia 


255.34 


s 


105.99 


32.23 


14.74 


100 


14 Irene 


151.30 


s 


105.39 


17.53 


5.98 


100 


48 Doris 


221.80 


c 


105.34 


1.58 


0.65 


98 


488 Kreusa 


150.12 


c 


101.99 


2.48 


1.18 


100 



" Density classes are assigned according to albedo and do not always correspond to taxonomies. 



MC experiments, we can compile a list of objects that should be 
modeled individually in an ephemeris. For N = 300 in the simul- 
taneous fit of section [53] a total of 523 asteroids were removed 
at least once from the test model during the 100 MC runs, and the 
average number of removed objects was approximately 240. In 
an ephemeris with an idealized asteroid model we should there- 
fore fit these 523 asteroid masses with the option of putting more 
than a half of the fitted values to zero. Among the 523 asteroids, 
there are 72 individuals removed from the test model on each 
run and 60 asteroids removed only once during the 100 runs. 
The distinction between asteroids having a high chance of re- 
moval and a low one is not clear, it is however possible to define 
an arbitrary limit above which the probability of being removed 



from the test model is reasonably high. Fixing this limit at 25% 
leads to a total of 287 objects, which are listed in Table Al of 
the supplementary online material. The table provides for each 
asteroid the probability of being selected as well as its maximum 
effects on the Earth-Mars, Earth- Venus, and Earth-Mercury dis- 
tances during the 1969 - 2010 time interval (for the standard 
mass set). It is interesting to note the existence of asteroids with 
relatively small effects on the Earth-Mars distance but with very 
high probabilities of being included in the individual part of the 
asteroid model: for example, 758 Mancunia induces a perturba- 
tion of approximately 10 m, but it is removed on each of the 100 
runs. A part of Table Al of the supplementary online material is 
reproduced in Table [5] 
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6.2. Systematic error estimation 

It is possible to use the results obtained in figure [TT to estimate 
systematic errors that will be induced by the residuals of the 
global perturbation during future missions like BepiColombo. 
These systematic errors can have a significant impact on the 
planned determinations of physical parameters from the ranging 
data. An extensive study of this problem is presented in Ashby 
|etaL| ( |2007> . 

The BepiColombo mission is expected to generate ranging 



model are its 99% modeling capacity and that, without the ring, 
systematic errors can reach several hundreds of meters. 



data to Mercury accurate down to 4.5 cm (see Ashby et al 
2007). Figure 1 1 shows that 4.5 cm per year is relatively close 



to the best possible residuals reached for the Earth-Mercury dis- 
tance. To take full advantage of this accuracy, the asteroid model 
used to process the ranging data will have to correctly account 
for approximately 200 individual asteroids. This is a relatively 
high number because obtaining accurate estimates of 200 aster- 
oid masses may still be difficult in the near future. If we estimate 
at 50 the number of asteroids that we are actually able to model 
with the highest accuracy, we can use figure 1 1 to obtain an es- 
timate of the systematic error for the BepiColombo mission if 
it took place today. For N = 50, fitting only the Earth-Mercury 
distance leads to residuals of approximately 6 m equivalent to a 
systematic error of 15 cm over a one year period. 



7. Discussion 

It was shown that with an appropriate selection scheme a ring is 
able to effectively reduce the amplitude of a perturbation induced 
by thousands of asteroids from an average 472 m to only 4 m. 
The ring thus represents more than 99% of the global perturba- 
tion, which clearly makes it a very suitable model. The value of 
99% was obtained with 100 MC experiments, so it can be con- 
sidered as quite robust. Also, this percentage does not depend 
on choices made for attributing asteroids with masses in the test 
model. The estimations of the amplitude of the residuals, the am- 
plitude of the global perturbation and the mass of the ring are, 
on the other hand, proportional to the mass of the test model and 
in consequence strongly dependent on the choices made in the 



attribution of asteroid masses. If for example in section 4.2 the 



interval used for the C density class were centered on 1 instead of 
1.5, all the previous parameters would have been approximately 
one third smaller. We are reasonably confident in the realism of 
the asteroid masses in our model because the mass of the ring ob- 
tain ed a fter removing the 300 most perturbing asteroids in sec- 

:0.15xl0- 10 M G 



tion 



5.1 



is relatively close to the value of 0.34 : 
obtained in INPOP06. 

The objective in this work was to show that the ring is a first- 
order model of a main-belt global effect and that it is able to rep- 
resent large numbers of objects in practice. The difficulty of fit- 
ting its mass with other highly correlated parameters is an impor- 
tant problem not considered here. In particular, the initial condi- 
tions of the planets were maintained fixed throughout our study, 
whereas they are fitted in an ephemeris. Because the global ef- 
fect acts mostly through linear drifts in mean-longitudes, a large 
part can be absorbed by changes of a few meters in the initial 
semi-major axes of the planets. The mass of the ring can be cor- 
related with other parameters as well, like the individual aster- 
oid masses or for example solar oblateness as shown by Fienga 
|et al.| l |2009) . These correlations can be considered as so impor- 
tant that the ring is eventually not implemented in the model (the 
case in DE421) or its mass is fixed to a certain value (the case 
for INPOP08). The major arguments for keeping the ring in the 



In section 2. 1 we fixed the radius of the ring to 2.8 AU and al- 
ways considered mass as the only parameter. We have briefly in- 
vestigated the possibility of fitting the radius and mass together. 
We find out that, when considering only the Earth-Mars distance, 
any change in mass can effectively be compensated for by a 
change in radius. In terms of the residuals on the Earth-Mars 
distance, moving the ring from 2.8 AU to 2.4 AU is equivalent 
to doubling the ring's mass. Similarly moving the ring to 3.4 AU 
is equivalent to dividing the mass of the ring by two. The same 
residuals can be obtained no matter the radius. Nevertheless in 
section [2~T| the mass of a ring with radius 2.8 AU indeed corre- 
sponds to the total mass of the represented asteroids (it is not the 
case for other radii). 

Throughout this paper, we estimated the amplitudes of per- 
turbations with the maximum reached on the 1969 - 2010 inter- 
val. This corresponds to the maximum norm || . Measuring 
perturbations in terms of the root mean square (equivalent to the 
norm || . H2) would divide all our amplitude estimations by ap- 
proximately three. This would lead to much more relaxed de- 
mands on the asteroid models. In particular an accuracy below 2 
m over 10 years would be reached in figure 11 for = 100 
instead of N - 300. Similarly, the final systematic error for 
BepiColombo with a correct N = 50 asteroid model would drop 
from 15 cm per year to 5 cm per year. We showed that the num- 
ber of asteroids that need to be accounted for individually is no 
more than 300, if in a more optimistic perspective all our ampli- 
tude estimations can be divided by a factor of three, the number 
of asteroids to account for individually could be as low as 100. 

This study was restricted to the main-belt perturbations. 
Other Solar System objects potentially have impacts on the 
ephemerides that should be estimated and possibly accounted 
for. We can mention trans -Neptunian objects already imple- 
mented in the EPM ephemeris or the Trojan asteroids whose ef- 
fect could be significant and which are certainly not accounted 
for by a ring. Also we only considered effects on the inner plan- 
ets because they provide the majority of accurate data today. 
The effect of asteroids on the outer planets can be non-negligible 
and should be considered in future studies; especially, effects on 
Jupiter can be significant because of the various resonances with 
the main belt. However, such studies for the outer planets will 
have an impact on data only when accurate observations of the 
outer planets are available over a sufficient time span. 



8. Conclusion 

A ring is an implementation of an averaged orbit, which is a very 
good model of the perturbation induced by the main-belt aster- 
oids. After removing less than 300 objects from the main-belt, 
the ring is able to account for more than 99% of the remaining 
perturbation on all the inner planets. Since the amplitude of the 
global effect can reach several hundreds of meters in terms of 
the Earth-Mars distance, it is advisable to keep the ring in an 
ephemeris model of the Solar System. 



Appendix A: 3D perturbing force of an asteroid ring 

We derive here the expression of the gravitational force exerted 
by a ring of center O, radius r', and mass m on a point P of mass 
M. The chosen reference frame is (O, i,j,k), with k orthogonal 
to the ring's plane and i the unit vector in the direction of the 
projection of r = OP on the ring's plane. We denote with r' 
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Fig. A.l. The reference frame 



we obtain, after straightforward computation, the expression of 
F in terms of the complete elliptical integral of first (7C) and 
second (£) kind : 

2GmM 

F = « i a& ^ r + « ! ~ a W(P) " £ 08)) r'i] . 

7ra(l -p)T /z 

This expression is valid for an internal or external body P with 
r + r' . Indeed in this case we always have < a < 1, hence 
also < /J < 1, so the elliptical integrals are well defined by 
their hypergeometric series, although one may use approxima- 
tion formulas for a fast evaluation. When the problem is planar 
(7 = 0), these expressions become equivalent to the expressions 
given in |Krasinsky~etaL] ( pO^j l or |Fienga et aT1 ( [2008l >. 

Acknowledgements. The first author wishes to acknowledge interesting discus- 
sions with James Hilton (US Naval Observatory). This work was done with the 
financial support of the CNES and the French Ministry of Education. 



the radius vector of a point on the ring. The longitude angle of 
r' with origin at i is denoted v' (see figure A. 1 1. The potential 
exerted by the ring on P is 



V(r) = - 



GmM 



2n 



i) 



dv' 



By defining / as the angle of r with the plane of the ring, the 
potential can be rewritten as 



V(r)-- 
With 



GmM 
2n 



(r 2 + r' 2 



dv' 



■ 2rr cos / cos v') 



1/2" 



2rr' cos / 

2 7T~ ' 
r + r 



we obtain 



V(r) = - 



GmM 1 



2ji 



dv' 



y r 2 + r /2 2n Jo (1 - o-cos v') 1/2 
The above expression can be rewritten as 
GmM 2 



V(r) 



Vf " 
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where r K(J3) is the complete elliptical integraQof the first kind, 
and (3 and T are defined as 



2a 



T = r ■ + r ■ + 2r'rcosI. 



1 + a 

The force exerted on P is F = -V r V. If 



8rV cos / 

V r (T) = 2r + 2r'i V t (J3) = r + 4r' 



2 , /2 

r + r 



1 The elliptical integral of first (7C) and second (£) kind are defined 
as fWhittaker & Watson||l9Z7] : 



f 

Jo 



J I - P sin 2 x 



dx ; £()8) : 



1 -ySsin xdx. 



P. Kuchynka et al.: A ring as a model of the main belt in planetary ephemerides, Online Material p 1 

Table Al. Asteroids selected for the individual part of the asteroid model when simultaneously fitting ranging data on all four inner planets. As 
described in section|6] only asteroids with probability greater than 25% of being removed from the test model during the 100 MC experiments are 
listed. For each asteroid this probabilty is given together with the asteroid's diameter, density class (in the standard mass set) and correponding 
perturbations on the Earth-Mars, Earth- Venus, and Earth-Mercury distances in terms of max \AD\ on the 1969-2010 time interval. Density classes 
are assigned according to albedo (see section |4~2| l and do not always correspond to taxonomies. 



Asteroid name 


Diameter [km] 


Density class 


Earth-Mars [m] 


Earth- Venus [m] 


Earth-Mercury [m] 


prob. [%] 


4 Vesta 


468.30 


S 


11455.34 


164.42 


323.36 


100 


1 Ceres 


848.40 


S 


7829.72 


2656.87 


935.08 


100 


2 Pallas 


498.06 


S 


6612.41 


977.16 


379.08 


100 


324 Bamberga 


229.44 


C 


1820.88 


17.04 


19.46 


100 


10 Hygiea 


407.12 


C 


806.70 


132.18 


61.95 


100 


6 Hebe 


185.18 


S 


691.70 


54.66 


21.15 


100 


532 Herculina 


222.38 


S 


660.82 


24.04 


6.53 


100 


19 Fortuna 


213.47 


C 


572.41 


18.53 


9.37 


100 


51 Nemausa 


147.86 


M 


504.89 


33.92 


12.06 


100 


7 Iris 


199.84 


S 


489.73 


37.20 


12.85 


100 


9 Metis 


152.70 


S 


364.56 


4.82 


1.87 


100 


20 Massalia 


145.50 


s 


285.84 


20.90 


7.31 


100 


139 Juewa 


156.60 


c 


281.71 


5.07 


2.78 


100 


31 Euphrosyne 


255.90 


c 


279.07 


18.07 


7.00 


100 


16 Psyche 


253.16 


s 


264.08 


38.66 


13.92 


100 


52 Europa 


302.50 


c 


252.60 


43.76 


15.93 


100 


24 Themis 


218.44 


c 


234.96 


2.07 


1.25 


100 


63 Ausonia 


103.14 


s 


234.40 


4.15 


1.44 


100 


747 Winchester 


171.72 


c 


222.34 


27.55 


13.06 


100 


60 Echo 


60.20 


s 


165.80 


0.15 


0.15 


100 


78 Diana 


120.60 


c 


157.41 


1.49 


0.41 


100 


354 Eleonora 


155.16 


s 


152.73 


0.37 


1.88 


100 


41 Daphne 


174.00 


c 


149.44 


6.11 


4.25 


100 


3 Juno 


233.92 


s 


145.29 


15.88 


5.58 


100 


654 Zelinda 


127.40 


c 


142.71 


21.79 


13.53 


100 


5 Astraea 


119.06 


s 


139.25 


10.29 


3.93 


100 


18 Melpomene 


140.56 


s 


130.84 


18.96 


14.11 


100 


128 Nemesis 


188.16 


c 


128.69 


17.32 


7.71 


100 


192 Nausikaa 


103.26 


s 


128.46 


0.67 


2.26 


100 


1 1 Parthenope 


153.34 


s 


115.67 


5.73 


3.86 


100 


106 Dione 


146.60 


M 


114.27 


4.99 


1.58 


100 


105 Artemis 


119.08 


c 


113.20 


1.33 


0.59 


100 


23 Thalia 


107.54 


s 


112.92 


1.05 


0.56 


100 


372 Palma 


188.62 


c 


112.84 


3.79 


1.58 


100 


29 Amphitrite 


212.22 


s 


106.71 


43.13 


13.78 


100 


419 Aurelia 


129.00 


c 


106.21 


3.20 


1.23 


100 


15 Eunomia 


255.34 


s 


105.99 


32.23 


14.74 


100 


14 Irene 


151.30 


s 


105.39 


17.53 


5.98 


100 


48 Doris 


221.80 


c 


105.34 


1.58 


0.65 


98 


488 Kreusa 


150.12 


c 


101.99 


2.48 


1.18 


100 


451 Patientia 


224.96 


c 


99.02 


4.18 


1.03 


99 


344 Desiderata 


132.28 


c 


96.77 


2.26 


1.94 


100 


109 Felicitas 


89.44 


c 


95.92 


1.79 


1.08 


98 


405 Thia 


124.90 


c 


95.73 


5.82 


3.21 


100 


511 Davida 


326.06 


c 


87.22 


12.60 


2.87 


100 


129 Antigone 


161.20 


M 


80.69 


26.62 


6.71 


93 


83 Beatrix 


81.38 


M 


78.82 


2.93 


1.68 


86 


88 Thisbe 


200.58 


c 


76.54 


14.11 


2.41 


99 


145 Adeona 


151.14 


c 


71.58 


2.53 


2.96 


87 
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Table Al. Continued. 



Asteroid name 


Diameter [km] 


Density class 


Earth-Mars [m] 


Earth- Venus [m] 


Earth-Mercury [m] 


prob. [%] 


187 Lamberta 


130.40 


C 


70.84 


9.56 


2.01 


94 


98 Ianthe 


104.46 


C 


70.43 


2.10 


0.83 


99 


13 Egeria 


207.64 


c 


68.44 


4.06 


2.41 


97 


27 Euterpe 


109.61 


s 


67.04 


2.09 


0.49 


100 


25 Phocaea 


75.12 


s 


65.91 


1.73 


0.32 


100 


240 Vanadis 


103.90 


c 


64.32 


1.94 


0.92 


100 


356 Liguria 


131.32 


c 


63.12 


1.22 


0.49 


100 


230 Athamantis 


108.98 


s 


63.11 


1.44 


0.35 


100 


22 Kalliope 


181.00 


s 


61.26 


6.96 


1.05 


100 


704 Interamnia 


316.62 


c 


61.16 


7.87 


6.50 


100 


505 Cava 


107.98 


c 


61.10 


1.65 


0.73 


100 


194 Prokne 


168.42 


c 


60.99 


4.43 


3.98 


99 


56 Melete 


113.24 


c 


59.77 


4.60 


1.62 


99 


268 Adorea 


139.88 


c 


58.53 


4.58 


1.53 


79 


53 Kalypso 


115.38 


c 


56.04 


0.50 


0.36 


99 


585 Bilkis 


58.10 


c 


55.47 


0.14 


0.20 


75 


94 Aurora 


204.88 


c 


54.04 


4.16 


3.04 


85 


74 Galatea 


118.70 


c 


53.67 


1.08 


0.78 


87 


259 Aletheia 


178.60 


c 


53.47 


1.90 


1.76 


100 


12 Victoria 


112.76 


s 


50.57 


12.31 


3.76 


100 


42 Isis 


100.20 


s 


50.15 


0.70 


0.50 


100 


141 Lumen 


131.04 


c 


49.65 


10.11 


5.24 


98 


8 Flora 


135.88 


s 


48.72 


3.24 


3.13 


100 


393 Lampetia 


96.90 


c 


46.80 


1.89 


1.54 


99 


81 Terpsichore 


119.08 


c 


46.67 


5.81 


3.29 


94 


87 Sylvia 


260.94 


c 


46.17 


15.28 


4.24 


100 


40 Harmonia 


107.62 


s 


45.73 


0.57 


1.49 


100 


694 Ekard 


90.78 


c 


43.86 


0.84 


0.74 


100 


37 Fides 


108.34 


s 


43.28 


0.31 


0.72 


100 


144 Vibilia 


142.38 


c 


43.26 


1.75 


1.05 


99 


36 Atalante 


105.60 


c 


42.90 


1.84 


0.83 


84 


410 Chloris 


123.56 


c 


41.81 


0.85 


0.24 


90 


313 Chaldaea 


96.34 


c 


39.34 


0.55 


0.92 


93 


50 Virginia 


99.82 


c 


39.16 


6.28 


2.36 


98 


173 Ino 


154.10 


c 


39.14 


2.81 


0.54 


100 


449 Hamburga 


85.60 


c 


39.06 


0.85 


0.15 


95 


516 Amherstia 


73.10 


s 


38.66 


0.53 


0.24 


99 


387 Aquitania 


100.52 


s 


38.34 


4.69 


2.14 


92 


164 Eva 


104.88 


c 


36.51 


1.57 


0.94 


98 


59 Elpis 


164.80 


c 


36.49 


8.28 


1.84 


72 


107 Camilla 


222.62 


c 


36.35 


13.62 


4.20 


93 


30 Urania 


100.16 


s 


35.99 


1.26 


0.14 


100 


247 Eukrate 


134.44 


c 


35.46 


0.57 


0.73 


66 


409 Aspasia 


161.62 


c 


35.33 


6.27 


2.04 


93 


17 Thetis 


90.04 


s 


35.24 


4.19 


1.59 


81 


386 Siegena 


165.00 


c 


34.72 


2.05 


0.75 


100 


49 Pales 


149.80 


c 


34.10 


9.27 


3.36 


97 


423 Diotima 


208.78 


c 


33.98 


10.06 


3.21 


100 


97 Klotho 


82.82 


s 


33.84 


0.35 


0.44 


94 


404 Arsinoe 


97.70 


c 


33.43 


0.99 


0.39 


66 


1021 Flammario 


99.40 


c 


32.21 


4.11 


3.02 


76 


46 Hestia 


124.14 


c 


31.98 


2.67 


1.65 


60 


135 Hertha 


79.24 


s 


31.96 


0.67 


1.10 


80 
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Table Al. Continued. 



Asteroid name 


Diameter [km] 


Density class 


Earth-Mars [m] 


Earth- Venus [m] 


Earth-Mercury [m] 


prob. [%] 


115 Thyra 


79.82 


S 


31.89 


0.66 


0.13 


100 


521 Brixia 


115.66 


C 


31.60 


6.43 


3.65 


100 


444 Gyptis 


163.08 


C 


31.31 


4.86 


2.37 


100 


212 Medea 


136.12 


c 


30.92 


4.67 


1.36 


85 


85 Io 


154.78 


c 


30.56 


5.84 


4.05 


99 


304 Olga 


67.86 


c 


30.50 


0.17 


0.14 


38 


21 Lutetia 


95.76 


s 


30.47 


2.31 


0.68 


98 


181 Eucharis 


106.66 


s 


30.32 


2.54 


0.86 


99 


89 Julia 


151.46 


s 


30.23 


17.04 


9.29 


77 


111 Ate 


134.56 


c 


29.28 


0.53 


0.18 


91 


471 Papagena 


134.18 


s 


28.89 


4.00 


0.90 


85 


96 Aegle 


170.02 


c 


28.86 


3.46 


2.34 


99 


44 Nysa 


70.64 


s 


28.79 


2.58 


0.94 


99 


416 Vaticana 


85.46 


s 


28.74 


0.96 


0.31 


79 


335 Roberta 


89.06 


c 


28.58 


0.96 


0.31 


83 


150Nuwa 


151.12 


c 


28.09 


2.91 


0.26 


84 


275 Sapientia 


96.68 


s 


27.62 


0.87 


0.11 


52 


216 Kleopatra 


135.06 


s 


27.33 


13.77 


4.38 


100 


69 Hesperia 


138.12 


s 


26.99 


5.74 


2.24 


97 


469 Argentina 


125.56 


c 


26.36 


1.38 


0.66 


48 


365 Corduba 


105.92 


c 


26.06 


3.75 


1.35 


92 


593 Titania 


75.32 


c 


26.01 


0.10 


0.12 


76 


680 Genoveva 


83.92 


c 


25.67 


0.86 


0.31 


99 


554 Peraga 


95.88 


c 


25.41 


7.98 


3.31 


91 


75 Eurydike 


55.90 


s 


25.22 


0.23 


0.13 


97 


130Elektra 


182.24 


c 


25.12 


4.74 


2.04 


97 


38 Leda 


115.94 


c 


23.89 


2.44 


1.07 


99 


200 Dynamene 


128.36 


c 


23.33 


3.63 


1.58 


100 


481 Emita 


108.47 


c 


23.02 


0.69 


0.22 


98 


914Palisana 


76.60 


M 


22.97 


3.39 


0.89 


94 


760 Massinga 


71.30 


s 


22.48 


0.93 


0.29 


94 


163 Erigone 


72.62 


c 


22.43 


1.31 


0.64 


95 


385 Ilmatar 


91.54 


s 


22.36 


3.52 


1.13 


51 


93 Minerva 


141.56 


c 


22.28 


1.96 


0.45 


53 


788 Hohensteina 


103.68 


c 


21.88 


0.56 


0.38 


100 


776 Berbericia 


151.16 


c 


21.86 


2.46 


0.53 


62 


702 Alauda 


194.72 


c 


21.38 


8.03 


3.15 


95 


193 Ambrosia 


65.97 


s 


21.31 


1.31 


1.00 


37 


455 Bruchsalia 


84.40 


c 


20.54 


1.27 


0.54 


97 


211 Isolda 


143.18 


c 


20.41 


1.41 


0.17 


90 


238 Hypatia 


148.50 


c 


19.91 


6.49 


1.94 


84 


762 Pulcova 


137.08 


c 


19.90 


0.59 


0.28 


100 


132 Aethra 


42.86 


s 


19.75 


0.48 


0.29 


83 


814Tauris 


109.56 


c 


19.71 


0.78 


1.10 


82 


415 Palatia 


76.34 


c 


19.51 


0.62 


0.57 


97 


626 Notburga 


100.74 


c 


19.20 


4.48 


1.87 


78 


674 Rachele 


97.34 


s 


18.85 


2.96 


1.12 


64 


28 Bellona 


120.90 


s 


18.81 


1.27 


0.68 


100 


68 Leto 


122.56 


s 


18.51 


3.28 


1.59 


99 


70 Panopaea 


122.18 


c 


18.38 


1.61 


0.76 


97 


751 Faina 


110.50 


c 


18.38 


1.76 


0.39 


68 


752 Sulamitis 


62.78 


c 


18.38 


0.48 


0.21 


26 


287 Nephthys 


67.60 


s 


18.22 


1.73 


0.61 


95 
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Table Al. Continued. 



Asteroid name 


Diameter [km] 


Density class 


Earth-Mars [m] 


Earth- Venus [m] 


Earth-Mercury [m] 


prob. [%] 


91 Aegina 


109.82 


C 


17.97 


3.53 


1.70 


63 


346 Hermentaria 


106.52 


s 


17.76 


4.07 


1.98 


84 


895 Helio 


141.90 


c 


17.48 


1.22 


0.34 


97 


322 Phaeo 


70.84 


c 


16.98 


0.28 


0.14 


46 


779 Nina 


76.62 


s 


16.92 


0.26 


0.45 


89 


350 Ornamenta 


118.34 


c 


16.55 


0.47 


0.18 


98 


34 Circe 


113.54 


c 


15.70 


0.14 


0.11 


49 


602 Marianna 


124.72 


c 


15.56 


0.24 


0.40 


94 


712 Boliviana 


127.56 


c 


14.91 


0.84 


1.41 


96 


325 Heidelberga 


75.72 


M 


14.75 


1.54 


0.50 


87 


786 Bredichina 


91.60 


c 


14.65 


0.69 


0.11 


34 


80 Sappho 


78.40 


s 


14.59 


3.96 


1.58 


76 


84 Klio 


79.16 


c 


14.36 


0.80 


0.29 


86 


308 Polyxo 


140.68 


c 


14.24 


0.47 


0.22 


60 


209 Dido 


159.94 


c 


14.02 


4.95 


2.13 


59 


241 Germania 


168.90 


c 


14.02 


2.63 


2.01 


87 


236 Honoria 


86.20 


s 


13.91 


0.89 


0.26 


73 


804 Hispania 


157.58 


c 


13.89 


5.51 


1.67 


75 


912Maritima 


83.18 


M 


13.81 


1.96 


1.05 


82 


337 Devosa 


59.10 


s 


13.75 


1.27 


0.38 


96 


54 Alexandra 


165.76 


c 


13.63 


9.30 


3.57 


93 


336 Lacadiera 


69.32 


c 


13.58 


0.84 


0.21 


90 


784 Pickeringia 


89.42 


c 


13.49 


0.30 


0.04 


93 


121 Hermione 


209.00 


c 


13.48 


3.85 


2.49 


67 


118Peitho 


41.72 


s 


13.25 


0.19 


0.09 


82 


266 Aline 


109.10 


c 


12.98 


0.08 


0.42 


54 


517 Edith 


91.12 


c 


12.24 


0.17 


0.07 


89 


72 Feronia 


85.90 


c 


12.21 


1.79 


0.56 


38 


595 Polyxena 


109.06 


M 


12.00 


2.94 


1.01 


37 


120 Lachesis 


174.10 


c 


11.81 


4.60 


0.30 


71 


599 Luisa 


64.88 


s 


11.80 


1.97 


0.95 


74 


112 Iphigenia 


72.18 


c 


11.71 


1.31 


0.73 


48 


61 Danae 


82.04 


s 


11.51 


1.20 


0.41 


80 


442 Eichsfeldia 


66.74 


c 


11.42 


1.49 


0.48 


50 


485 Genua 


63.88 


s 


11.30 


0.29 


0.14 


87 


45 Eugenia 


214.62 


c 


11.04 


11.47 


4.07 


100 


582 Olympia 


43.40 


s 


11.04 


0.17 


0.08 


48 


39 Laetitia 


149.52 


s 


10.75 


1.65 


2.26 


98 


568 Cheruskia 


86.98 


c 


10.74 


2.56 


0.99 


62 


156 Xanthippe 


120.98 


c 


10.70 


5.90 


1.03 


83 


148 Gallia 


97.76 


s 


10.53 


0.27 


0.17 


93 


1171 Rusthawelia 


70.12 


c 


10.43 


0.39 


0.17 


47 


134 Sophrosyne 


123.26 


c 


10.30 


1.70 


0.27 


84 


758 Mancunia 


85.48 


s 


10.30 


0.48 


0.17 


100 


345 Tercidina 


94.12 


c 


10.25 


0.58 


0.12 


73 


667 Denise 


81.28 


c 


10.13 


0.09 


0.05 


57 


323 Brucia 


35.82 


s 


10.10 


0.15 


0.07 


63 


92 Undina 


126.42 


s 


9.91 


3.70 


1.96 


44 


79 Eurynome 


66.46 


s 


9.86 


0.44 


0.17 


41 


221 Eos 


103.88 


s 


9.79 


1.60 


0.51 


93 


283 Emma 


148.06 


c 


9.79 


5.34 


1.81 


81 


503 Evelyn 


81.68 


c 


9.67 


0.09 


0.05 


50 


618 Elfriede 


120.30 


c 


9.63 


0.98 


0.42 
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Table Al. Continued. 



Asteroid name 


Diameter [km] 


Density class 


Earth-Mars [m] 


Earth- Venus [m] 


Earth-Mercury [m] 


prob. [%] 


47 Aglaja 


126.96 


C 


9.54 


1.96 


1.46 


66 


172 Baucis 


62.42 


s 


9.43 


0.09 


0.14 


95 


77 Frigga 


69.24 


s 


9.37 


0.75 


0.36 


97 


566 Stereoskopia 


168.16 


c 


9.28 


1.37 


0.47 


80 


217 Eudora 


66.24 


c 


9.24 


0.50 


0.27 


77 


171 Ophelia 


116.68 


c 


9.06 


3.65 


1.21 


71 


739 Mandeville 


107.54 


c 


8.94 


0.39 


0.46 


69 


790 Pretoria 


170.38 


c 


8.83 


0.43 


0.37 


51 


147 Protogeneia 


132.94 


c 


8.64 


3.24 


0.95 


58 


86 Semele 


120.56 


c 


8.25 


1.28 


0.78 


30 


62 Erato 


95.40 


c 


8.24 


0.75 


0.23 


70 


250 Bettina 


79.76 


s 


8.24 


1.12 


0.41 


56 


71 Niobe 


83.42 


s 


8.09 


0.26 


0.25 


49 


185 Eunike 


157.50 


c 


7.88 


5.32 


2.08 


43 


591 Irmgard 


51.86 


c 


7.80 


0.10 


0.06 


36 


234 Barbara 


43.76 


s 


7.69 


0.16 


0.06 


93 


978 Aidamina 


78.74 


c 


7.68 


0.16 


0.06 


47 


679 Pax 


51.46 


s 


7.64 


1.02 


0.34 


83 


420 Bertholda 


141.24 


c 


7.39 


3.05 


1.02 


57 


95 Arethusa 


136.04 


c 


7.35 


0.72 


0.50 


63 


735 Marghanna 


74.32 


c 


7.34 


2.80 


1.58 


63 


769 Tatjana 


106.44 


c 


7.23 


0.30 


0.14 


57 


210 Isabella 


86.66 


c 


7.15 


0.48 


0.11 


42 


117 Lomia 


148.72 


c 


7.09 


2.50 


0.30 


45 


491 Carina 


97.30 


c 


7.08 


1.23 


0.54 


64 


1963 Bezovec 


44.68 


c 


6.96 


0.09 


0.02 


62 


122 Gerda 


81.68 


s 


6.88 


0.59 


0.23 


77 


196 Philomela 


136.38 


s 


6.88 


5.20 


1.88 


54 


388 Charybdis 


114.18 


c 


6.84 


1.39 


0.11 


52 


223 Rosa 


87.60 


c 


6.73 


0.89 


0.50 


58 


1107 Lictoria 


79.18 


c 


6.73 


0.12 


0.06 


41 


431 Nephele 


95.04 


c 


6.69 


0.87 


0.38 


68 


168 Sibylla 


148.38 


c 


6.65 


1.04 


0.27 


41 


1036 Ganymed 


31.66 


s 


6.49 


6.92 


3.74 


76 


138 Tolosa 


45.50 


s 


6.43 


0.37 


0.14 


73 


949 Hel 


69.18 


c 


6.34 


0.33 


0.11 


35 


893 Leopoldina 


76.14 


c 


6.29 


0.70 


0.26 


49 


179 Klytaemnestra 


77.68 


s 


6.27 


0.88 


0.34 


78 


675 Ludmilla 


72.09 


s 


6.24 


0.38 


0.08 


76 


445 Edna 


87.18 


c 


6.18 


0.17 


0.19 


38 


718 Erida 


72.94 


c 


6.13 


0.79 


0.23 


32 


162 Laurentia 


99.10 


c 


6.10 


1.03 


0.62 


35 


375 Ursula 


182.53 


c 


6.05 


5.83 


1.50 


64 


381 Myrrha 


120.58 


c 


5.86 


0.17 


0.14 


49 


1001 Gaussia 


74.68 


c 


5.68 


0.24 


0.12 


36 


306 Unitas 


46.70 


s 


5.64 


0.43 


0.15 


72 


690 Wratislavia 


134.66 


c 


5.63 


4.23 


1.38 


61 


1013 Tombecka 


31.92 


s 


5.62 


0.01 


0.02 


26 


1694 Kaiser 


29.06 


c 


5.62 


0.16 


0.07 


35 


76 Freia 


183.66 


c 


5.61 


2.84 


1.24 


44 


55 Pandora 


66.70 


s 


5.57 


0.46 


0.04 


76 


43 Ariadne 


65.88 


s 


5.43 


1.97 


0.91 


44 


604 Tekmessa 


65.16 


c 


5.43 


0.45 


0.30 
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Table Al. Continued. 



Asteroid name 


Diameter [km] 


Density class 


Earth-Mars [m] 


Earth- Venus [m] 


Earth-Mercury [m] 


prob. [%] 


253 Mathilde 


58.04 


C 


5.40 


0.36 


0.24 


41 


202 Chryseis 


86.16 


s 


5.33 


1.58 


0.58 


39 


569 Misa 


72.96 


c 


5.31 


0.28 


0.08 


32 


176 Iduna 


121.04 


c 


5.20 


0.69 


0.52 


49 


26 Proserpina 


94.80 


s 


5.08 


3.28 


1.08 


52 


33 Polyhymnia 


53.68 


s 


5.03 


0.20 


0.23 


38 


67 Asia 


58.10 


s 


4.98 


0.73 


0.29 


59 


791 Ani 


103.52 


c 


4.97 


0.38 


0.06 


38 


980 Anacostia 


86.18 


s 


4.83 


1.80 


0.61 


54 


165 Loreley 


154.78 


c 


4.76 


0.77 


0.54 


31 


1015 Christa 


96.94 


c 


4.72 


0.41 


0.28 


26 


971 Alsatia 


63.76 


c 


4.71 


0.49 


0.22 


30 


583 Klotilde 


81.64 


c 


4.65 


0.30 


0.02 


28 


952 Caia 


81.60 


c 


4.64 


1.10 


0.29 


31 


663 Gerlinde 


100.88 


c 


4.51 


0.31 


0.07 


26 


584 Semiramis 


54.00 


s 


4.47 


0.32 


0.30 


42 


152 Atala 


122.84 


s 


4.36 


2.63 


0.95 


27 


349 Dembowska 


139.78 


s 


4.00 


3.56 


1.87 


41 


127 Johanna 


124.55 


s 


3.98 


0.69 


0.31 


26 


178 Belisana 


35.82 


s 


3.93 


0.11 


0.01 


25 


204 Kallisto 


48.58 


s 


3.91 


0.49 


0.15 


29 


65 Cybele 


237.26 


c 


3.86 


1.83 


0.57 


65 


126 Velleda 


44.82 


s 


3.75 


0.30 


0.12 


26 


328 Gudrun 


122.92 


c 


3.75 


1.63 


1.15 


42 


57 Mnemosyne 


112.60 


s 


3.60 


3.22 


1.65 


79 


665 Sabine 


51.10 


s 


2.89 


0.02 


0.03 


25 



